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Abstract 

In this paper the mode III crack problem for two bonded homogeneous half planes is 
considered. The interfacial zone is modelled by a nonhomogeneous strip in such a way that 
the shear modulus is a continuous function throughout the composite medium and has 
discontinuous derivatives along the boundaries of the interfacial zone. The problem is 
formulated for cracks perpendicular to the nominal interface and is solved for various crack 
locations in and around the interfacial region. The asymptotic stress field near the tip of a 
crack terminating at an interface is examined and it is shown that, unlike the corresponding 
stress field in piecewise homogeneous materials, in this case the stresses have the standard 
square root singularity and their angular variation is identical to that of a crack in a 
homogeneous medium. With application to the subcritical crack growth process in mind, the 
results given include mostly the stress intensity factors for some typical crack geometries and 
various material combinations. 

1. Introduction 

In studying the fracture process in bonded materials, there are two groups of problems 
in which the mechanical modelling of interfacial regions may play a major role. The first is 
”debonding” and the second, crack penetration. The investigation of the debonding process 
requires the solution of a crack problem in which the crack is located in the interfacial zone 
and is parallel to the nominal interface. To study the second problem one may need the 
solution of a crack which lies in a plane generally perpendicular to the nominal interface and 
which intersects the interfacial region. Even though there are cases such as certain geophysical 
problems and some transient thermal stress problems involving materials with heavily 
temperature dependent material properties in which the constituent materials under 
consideration are inherently nonhomogeneous, in many engineering applications the logical 
assumption has been to treat the bonded materials as piecewise homogeneous media. Within 
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the context of the linear elastic theory the only time such an idealized model of the interfacial 
regions leads to certain physical anomalies appears to be the presence of a crack in the medium 
that either lies along or intersects the interface. 

As observed, for example, by Muskhelishvili (1953), Erdogan (1963,1965), England 
(1965), and Rice and Sih (1965), in the case of an interface crack the power of the crack tip 
singularity is complex predicting a physically unacceptable result, namely the interference of 
crack surface displacements (see also Comninou (1977) and Atkinson (1977) for further 
discussion of the problem based on the crack closure concept). For a crack intersecting the 
ideal bimaterial interface perpendicularly the physical difficulty arises from the fact that the 
analysis gives a singularity that is not the standard square root type. It was shown, for 
example, by Cook and Erdogan (1972) that for a crack terminating at a bimaterial interface, 
depending on the relative magnitudes of material constants, the power of stress singularity 
may be greater or less than 0.5. Similarly, Erdogan and Biricikoglu (1973) showed that for a 
crack crossing the interface the stress state exhibits a weak power singularity at the point of 
intersection (see also Bogy (1968,1971) and Hein and Erdogan (1971) for results in bonded 
wedges). The consequence of these analytical results is that in such problems, since the 
stresses and displacements near the crack tip do not remain self-similar as the crack 
propagates, it is not possible to calculate the strain energy release rate and to apply the 
techniques of the conventional fracture mechanics. 

For a more accurate representation of the actual physical problem in bonded materials 
involving cracks it is clear that one needs to reexamine both the material model and the 
interface model. Maintaining the ideal interface model, Knowles and Sternberg (1975) and 
Shih and Asaro (1987) studied the interface crack problem for a nonlinear elastic material 
bonded to a rigid half space. Even though the anomalous deformation behavior at the crack 
tip may perhaps be rectified by a proper choice of the material model, it is highly unlikely that 
changes in material model alone will be sufficient to affect the standard square root singularity 
for a crack terminating at the interface and to eliminate the weak power singularity for a crack 
crossing the interface. This is then an important motivation to experiment with the 
mechanical modelling of interfacial zones. Another reason is, of course, the actual physical 
nature of interfacial zones. In many cases, regardless of the particular mechanism of the 
binding process, for thermodynamical consistency the interfacial region needs to be considered 
as a distinct phase (see, for example, the review article by Good (1967)). There is always 
some interdiffusion taking place across the interface so that there is usually a gradual change 
in the volume concentration of certain elements belonging to the adjacent material. Because of 
this even in diffusion bonding in which the thickness of the reaction zone may be only a few 
atomic layers, the actual thickness of the region having very steep variation in material 
properties may be substantially greater. 
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The basic assumption made in this study in modelling the interfacial zone, therefore, is 
that the mechanical properties of the composite medium may exhibit no jump discontinuties. 
The interfacial zone is thus aassumed to be a thin elastic layer between two homogeneous 
materials with steeply varying mechanical properties. The corresponding plane elasticity 
problem for a crack lying in the nonhomogeneous interfacial region and parallel to the nominal 
interface was considered by Delale and Erdogan (1988a). The important special case of the 
’’interface” crack in which the crack is located in the plane where the variation of the Young’s 
modulus forms a kink was studied by Delale and Erdogan (1988b) who showed that by using 
this model the oscillatory behavior predicted by the ideal interface model does indeed disappear 
(see also Atkinson (1977) and Atkinson and List (1978) for related studies). The main interest 
in this study is in the crack penetration problem. To cover all possible crack configurations it 
is assumed that the medium contains three collinear cracks one in each constitutent material 
oriented perpendicularly to the nominal interface. The medium is assumed to be under 
antiplane shear loading with displacements specified away from the crack region. The 
corresponding crack problem for two bonded nonhomogeneous half planes under mode III 
loading conditions was studied by Erdogan (1985). The important result found by Erdogan 
was that for the crack tip terminating at the kink plane of modulus distribution the stress 
state has the standard square root singularity. 

2. Formulation of the Problem 

Consider the antiplane shear problem shown in Fig. 1. Assume that the medium is 
loaded away from the crack region, the ’’homogeneous” problem in the absence of cracks has 
been solved, and 

fl’jyzC** 0 ) = T j( x ) = *Oyz(x,0) , aj<x<bj , i=l,2,3 (1) 

are the only external forces. We further assume that, compared to the crack lengths and the 
layer thickness h, the thickness of the homogeneous materials 2 and 3 are large and hence can 
be treated as being infinite, and the shear modulus of the interfacial layer is given by 

/i 1 (x) = /i 2 e^ x , 0<x<h . (2) 

From (2) and /^(h) = n 3 it follows that 

P - h lo g(/^ 3 /P2) • (3) 

Noting that y=0 is a plane of symmetry, the problem shown in Fig. 1 may be formulated as 
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follows: 


V 2 w x + 0 = 0 , 0<x<h , 0<y<oo , 

V 2 w 2 = 0 , -oo<x<0 , 0<y<oo , 

V 2 w 3 = 0 , h<x<oo , 0<y<oo , 

w x (0,y) = w 2 (°,y) , ff lx2 (0,y) = <r 2xz (0,y) , 0<y<oo , 
w i(h*y) = w 3 (h,y) , o-i xz (h,y) = 0 - 3xz (h,y) , 0<y<oo , 
w 1 (x,0) = 0, 0<x<aj, b 1 <x<h , 

^ly 2 ( x »°) = r l( x ) » a 1 <x<b 1 , 

w 2 (x,0) = 0 , -oo<x<a 2 , b 2 <x<0 , 

<r 2yz( x »°) = t 2( x ) » a 2 <x<b 2 , 

w 3 (x,0) = 0 , h<x<a 3 , b 3 <x<oo 


^yzC** 0 ) = r 3( x ) * a 3 <x<b 3 . 

A 

We now express the solution of (4)-(5) in the following form : 


r 00 f 00 

vj( x ,y) = fj(y,a)e*'" x da + f 9 j(x,a)sin(ya)da + Cj, (j=l,2,3). 

» -oo J O 


From (4)-(6) and (12) it may be shown that 


(4) 

(5) 

( 6 ) 
(7a, b) 
(8a, b) 

(9a, b) 
(10a, b) 
(11a, b) 
( 12 ) 


2 In (12) Cj_, c 2 , c 3 are included to account for possible rigid body displacements when 
the cracks cross the interfaces. If, however, in the analysis the continuity conditions (7a) and 
(8a) are replaced by 

J- w x (0,y) = ^ w 2 (0,y) , ^ w x (h,y) = w 3 (h,y) , 0<h<oo, (A) 

these constants disappear and conditions (9a), (10a) and (11a) along with (A) insures the 

single-valuedness and continuity of the displacements. 
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fl(y»«) = A 1 (a)e" liy , y>0 , (13) 

g x (x,a) = B 1 (a)e DlX + B 2 (a)e n2X , 0<x<h , (14) 

f 2(y> a ) = A 2 (a)e ^ y , y>0 , (15) 

g 2 (x,a) = D(a)e ax , x<0 , (16) 

f 3 (y> a ) - A 3 (a)e ^ V , y>0 , (17) 

g 3 (x,a) = C(a)e‘ ax , x>0 (18) 


where 


mi = -|oNl+^ , 


/?, P 

n l = - o+ Q l> n 2 = ■ o' 0 !’ °1 = 






Introducing the unknown functions 
<£j( x ) = ^ wj(x,0) , aj<x<bj , j= 1,2,3 


from (12)-(19) and (21) we obtain 
• r b j 

A j( Q ) = 5 j a <^j(t)e lat dt . 


(19) 

( 20 ) 


( 21 ) 


( 22 ) 


It should be observed that at y=^0 the solution must satisfy the following conditions (Fig. 1): 


lim 

y-+0 


w 2 (x,y) = 0 , x<a 2 , 


ylim +o w 2 (0,y) = y lmn +o w x (0,y) , 


lim , n w i( h »y) = lim , n w 3( h »y) » 

y->+0 y->+0 


lim w 3 (x,y) = 0 , x>b 3 . 

y-’+O 


(23a-d) 


By substituting from (12)-(18) and (22) into (23a-c) we find 
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1 f 2 

c 2 = 2 <t> 2 (t)dt , 

j a 2 


c l" c 2 = 2 } a ^l( t ) dt + 2 ^2( t ) dt * 

r bi *b 3 

c 3* c l = 2 + 2 > 

J a^ J a 3 


(24a-c) 


which determine the constants c 1? c 2 and c 3 . The fourth condition in (23) gives 


c 3 — "2 [ ^3 00 ^ 

J a 3 

which, with (24), implies that 

3 r b j 

£ ^(dt = 0 . 

1 Jaj J 


(25) 


(26) 


Note that if the cracks are disconnected (9a), (10a) and (11a) require <f > j defined by (21) to 
satisfy the following conditions: 

f b j 

*j(t)dt = 0, j=l,2,3. (27) 

J aj 


From (24) and (27) it is seen that in this case the constants Cj would be zero. However, if two 
or all three cracks are connected, then the appropriate combination of (26) and (27) provide 
the single- valuedness conditions and (24) gives c x , c 2 and c 3 . 

If we now substitute from (12)-(18) into the homogeneous conditions (7) and (8), 
observe that 


3wj 

^jxy = ^i~dx ’ 


0Wj 

'jy* = 


j=l,2,3 , 


(28) 


invert the sine transforms, we obtain 


B^a) + B 2 (a) - D(a) = R x (a) , 

n l®l( a ) "I" n 2®2( a ) ‘ a D( a ) = R 2 (a) , 
e Illll B 1 (a) + e^^B^a) - e" ah C(a) = R 3 (a) , 
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(29a-d) 


n l e 1 B 1 (a)+n 2 e 2 B 2 (a)+Qfe' a ^C(a) = R 4 (a) , 

where the functions Rj, (i=l,....,4) are given in Appendix A. Thus, from (29) and (A6) it 
follows that 

B l( .) = + <^2_V"> h e"^ h - , ) Wl («)d, 

*(a+n 2 )[ e - (a-n 2 ) [ e 1 , 

J a 2 * a 3 


B 


,<«) = i{2^ ) j^^V ni ‘ • <^e’”' h «” 2<h -* ) ]d > (.)dt 


+(a-(-n 1 )f e at ^ 2 (t)dt + (“*n 1 )[ e 1 e ^ 3 (t)dt) , 
J a 2 * 


3 -nih -a(t-h) 
e A “ 

3 


C(«) = S ?{|" 1 [^ ) « n2h '" ll ‘ • ^e Mh - ,) ]d 1 (.)d. 


+ 2a, f b2 «“» h e “‘d 2 (t)dt - I [ b2 (l. e - 2 “ lh )e a(t - | ‘ ) d3(‘)dt} , 

J a 2 * a*^ 


D(«) = l{f bl [.2teHl)e" ,l ‘ + 

J 


- | f b2 (l- e ' 2Qlh )e°V 2 (t)dt - 2a x fV “^e'^WtJdt} , (30a-d) 

z J a 2 Ja 3 


where 


A = (a-n 1 )(a+n 2 )e 2 “ lh - (a+n 1 )(a-n 2 ) 


(31) 


3. The Integral Equations 

In the previous section the problem has been formulated in terms of three unknown 
functions <f> v <j> 2 and <j> 3 . These functions are determined by using (9b), (10b) and (lib) which 
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are the only conditions that have not yet been satisfied. Thus, by substituting (12)-(18) and 
(28) into (9b), (10b) and (lib) we obtain 


yHm+o *iy Z (*>y) = ^ y lim +() | ^m 1 A 1 (a) e miy '' aX da 




J o 


X / X n 2 X n , /X 

+B 2 (a)e * ]ada=r x (x), a 1 <x<b 1 , 


yl^Vo <T 2yz(x»y) 


t2 

2?r 


lim 

y-*+o 


r 

V -c 


-|«|A 2 (a)e 


-|a|y-iax 


da 


+ 




da = r 2 (x), a 2 <x<b 2 , 


y S“ + 0 ff 3y 2 (x,y) 


& li m 

2* y _+o 


r oo 

I -MA 3 (a)( 

j -oo 


-|a|y-iax 


da 




da = r 3 (x) , a 3 <x<b 3 


(32a-c) 


By substituting from (22) and (30) into (32) we obtain the following system of integral 
equations for <f> j, (j=l ,2,3): 


lim 


y^+o 



by(x,y,t)^j(t)dt = Tj(x) , aj<x<bj . 


(33) 


The singular parts of the kernels in (33) can be separated after performing the 
appropriate asymptotic analysis. Note that the leading terms in the asymptotic expansion of 
the infinite intergrals giving the kernels in (33) are identical to that given by Erdogan (1985). 
Therefore, the details of the analysis will not be given in this paper. Thus, by following the 
procedure described by Erdogan (1985) (33) may be reduced to 



a 

t-x 


+ k jj(x,t)]<£j(t)dt = ^Tj(x) , aj<x<bj 


(34) 


where the kernels ky, (i j = l,2,3) are given in Appendix B. 

In examining the asymptotic behavior of the kernels for cracks fully embedded in the 
three materials shown in Fig. 1, that is for b 2 <0, a 1 >0, b x <h and a 3 >h, it is sufficient to 
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separate the singular parts of the kernels that become unbounded for x=t only. However, if 
one or two crack tips touch one of the interfaces x=0 or x=h, then it is necessary to separate, 
in addition to the parts that are singular at x=t, the terms that become unbounded at x=t=0, 
(a 2 <(t,x)<b 1 ) and at x=t=h, (a 1 <(t,x)<b 3 ). The inspection of the kernels given in (34) and 
Appendix B shows that for embedded cracks as expected the only unbounded terms in the 
kernels that may influence the nature of the crack tip singularity are (t-x)“^. The remaining 
unbounded terms are logarithmic (i.e., log|t-x|) which are square integrate and may be 
treated as Fredholm kernels. It is also seen that the terms that become unbounded for x=t=0 
and x=t=h are log|t-fx| and log|2h-x-t|. Again, these terms would have no influence on the 
crack tip singularity. One may, therefore, conclude that in this problem for all crack 
geometries the displacement derivative at the crack tip would have the standard square root 
singularity. Also, as in the case of piecewise homogeneous materials considered by Erdogan 
and Cook (1974), there would be no singularity at the point of intersection of the crack and 
the interface. 


4. Stress Distribution Near the Crack Tips and Stress Intensity Factors 

The stress distribution near the crack tips may be obtained by examining the 
asymptotic behavior of the corresponding expressions found by substituting from (12)-(18) and 
(30) into (28) for small values of r, where r is the radial distance from the crack tips. 
Consider, for example, the distribution of < 7 lyz around an embedded crack in material 1. From 
(32a) and (22) the leading term for ^ lyz may be expressed as 


r lyz( x »y) - ^ 2 ? [ ti(t)dt f 
J a^ J “C 


lm 

~a 


i mjy ia(t-x) 

- e e da 


( 35 ) 


Also, from (34) it is clear that the solution of the singular integral equations is of the form 


*j(t) = 


J(t-a-)(bj-t) 


aj<t.bj . 


(36) 


Thus, observing that for |a|—>oo m^-lal and defining 


_ bj-a^ ^ b 1 +a 1 * _ b 1 -a 1 ( ^ b 1 +a 1 ^ _ b 1 -a 1 


1 = - V r + “I 


i + 


y = 


a = b^i S) G i( t ) = F i(0 


(37) 


(35) becomes 
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( 38 ) 


r ( x y)— ItlW f°° 111 e' |s|,? ' is ^d 8 f 1 Fl(r) e isr dr 

r lyz( x »y)- ir(bl .a 1 )J. oo S « % ^2 • 


Again, since the possible unbounded terms in (38) are determined by the leading terms 
integrand for s—kx>, by using (see, for example, Gradshteyn and Eyzhik, (1965)) 


F j(r) 


,isr 


UP 


dr = 


N 2|s| 


-i( 8 -^) 

[F x (l)e ( 48 ) +F 1 (-l)e 4s 1+0(1) 


we find 


_ , x * 2 ^i( x ) ri r 

^iy 2 (x,y) - jr (b i . ai ) \| 2 J Q 


«> -ns 

S^-{-F 1 (l)sin[((-l)s+|] - Fj^lJsin^^+lis-^Jds. 


Similarly from (see Gradshteyn and Ryzhik, (1965)) 


f a ^ ^ ( "1 I> = ^^73 { £jo-- 1 » 


and (40) it may easily be shown that 


lyz 


(x y) ^ nM 10 

{ y> ~ (b r a0/2 $T P 


cos | + 


A*i(*) Fi(-i) e 

(b^F2 -Jzf Sin 2 


where 


pco&B = £-1, psin0 = t) and />cos0 = £+1, psin0 = r ? 


near £=1 and £=-l, respectively. If we now use (37), from (42) we obtain 

*lyz( x >y) = ~ j ( b^ai)/2 ^ TCOSe cos I , rcos<? = x-b x , rsin* = y 

near x=bj, y=0 and 


0- lyz( x »y) - I e /? rcos 0 sin ^ Icoa0 _ x . ai , rsinfl = y , 

^(b 1 - a i)/2 


near x=a x , y=0. 

Based on the asymptotic results expressed in (44) and (45) one may observe that 


in the 

(39) 

(40) 

(41) 

(42) 

(43) 

(44) 

(45) 
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(a) If we define the mode III stress intensity factors 


M a i)=l‘m ai ^ 2(a x -x) <r lyz (x,0), ksCb^sliin^ J2(x-b x ) * lyz (x,0) , (46) 


from (44) and (45) it is seen that 


Mll) = °>w • Mb,) = G l( b,) . 


(47a, b) 


These results are identical to that found for homogeneous materials* 

(b) Note that since expression (32a) for a lyz is valid for 0<x<h and since the form of 
the solution given by (36) for <£ x (t) remains unchanged for a x =0 and b x =h, the expressions 
giving the stress distribution around the crack tips would remain valid for the limiting cases 
a x =0 and bj=h provided 0 is restricted to 0<6<n/2 for a x =0 in (45) and ir/2<0<n for b x =h 
in (44). 

(c) Expressions similar to (44) and (45) can easily be developed for ^i xz (x,y) near the 
crack tips. It can indeed be shown that the expressions for <r lxz are identical to (44) and (45) 
provided cos(0/2) is replaced by sin(0/2) in (44) and sin(0/2) by cos(0/2) in (45). 

(d) Aside from the factor e ^ rcos ^ } the stress distributions (44) and (45) are identical to 
that found for homogeneous materials. However, note that the expressions are valid only for 
“small” values of r for which e ^ rcos ^ ; s ne arly unity. 

(e) From the foregoing results it does not automatically follow that for the limiting 
cases a x =0 and b x =h (Fig. 1) (47) is still valid and the validity of the expressions (44) and 
(45) can be extended into materials 3 and 2, respectively. To examine this problem we now 
consider the distribution of <r 2 yz (x,y) near the crack tip (a x , 0 ) for the limiting case of a x =0 
b 2 <0. Following Erdogan (1985), the leading term for <r 2yz may be expressed as 3 


r 2 yzV 


/ \ f ^ 1 , / x (/?t-7x)/2 f 00 -a(t-x) 

K (x,y) = -2^ — <£ x (t)e v } dt e v } cosay da . 

Jo Jo 

By using the definitions (37) (with F 1 (r)=:G 1 (t)e^ t ^) from (48) we find 

, \ , 2/i 2 e Tx/2 f°° s£ , T 1 F i(r) 

*2yz(*#) = e cosrjs ds e dr. 

1 J -1 ^ 1-T 


(48) 


(49) 


yx Qx 

3 In Erdogan (1985) /i 2 as well as /i x is a function of x (fi 2 e anc ^ PlF > respectively). 
Thus, by assuming that both f3 and 7 are nonzero, the asymptotic results given here cover all 
limiting crack geometries that may arise from Fig. 1. 
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If we define 


oo 


F l( r ) = E A n T n (r) , (T n (r) = cosn0 o , T = cos 6 0 ) 

0 


(50) 


and use (see Erdelyi (1953)) 

f e scos ^o cogn ^ o( j^ 0 -s (-i) n x I n (s) 

Jo 

(49) becomes 


(51) 


vx/2 


a 2yz( x >y) - £(-0 nA n [ e S ^I n (s)cos»jsds 

0 ■'o 


(52) 


where T n is the Chebyshev polynomial and I n is the modified Bessel function. Now by 
observing that 


s oo 

Tn (s) = q== for s->oo and £ (-J) A„ = Fj(-l) , 


(53) 


from (37), (52) and (53) we find 

ff2yz( x >y) = eTX ^ F l(‘ 1 ){ o e ax cosayda 


(54) 


In material 2 near the crack tip (a 1 =0, y=0) we have 


x<0 , -x = rcosflj , y=rsin^ 1 , 


(55) 


where 6^ is measured from the negative x axis. Thus from (54), (41), (55) and F 1 (-1)=G 1 (0) 
it may be shown that 

-^7rcos $ 1 ^ i/?rcos^ 1 


T /„ V N - A*2 e " G 1 (rcos0 1 )e' 

< ' 2 » <J! ’ y) - ^ TIT 


e , 


cos 


2 • 


(56) 


If we define the stress intensity factor at the crack tip (a 1 =0, y=0) by 


M°) 0 <r 2yz( x, °) 


(57) 
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from (56) and (55) it is seen that 


MO) 



Gx(0) , 


( 58 ) 


which is identical to (47a) with a x =0. Since 0 X =7 r-0, also note that aside from the factors 
exp(/?rcos0), exp(^rcos0 1 ) and exp(-^7rcos0 1 ) the stress distributions (45) and (56) are 
identical. From a practical viewpoint since all asymptotic expressions such as (56) are valid 
only for r— >0, in (44), (45) and (56) all factors involving e r should be ignored. This leads to a 
very simple and important conclusion, namely that the behavior of the singular field near the 
crack tip in nonhomogeneous materials is identical to that in homogeneous materials provided 
/i(x,y) is a continuous (but not necessarily a differentiable) function. 


5. The Solution and Results 

By defining the normalized quantities 


t = 


bj-a : 

T 


r + 


bj+ar 


x = 


V a j 


s + 


bj+aj 


^(t) = fj(r) 


V a j 


2 k ij( x ’ t ) = K u( r ’ s ) > ^7(7) r i( x ) = 8i( s )» (ij=1.2,3), 


(59) 


for disconnected cracks the integral equations (34) and conditions (27) may be expresed as 
rl 3 6 ;; 


I E [f4 + Kjj(s,r)]fj(r)dr = gj(s) , -1<x<1, (i=l,2,3) , 

J -1 1 


(60) 


| fj(r)dr = 0 , (j=l,2,3) . 
Equations (60) are solved by letting 



(61) 


(62) 


and by using 


Tn< ‘"f <i - C) = a-Un-lC 8 ) , (n>l) , |s|<l , 
-1 (r-s ) n 1-r 2 


(63) 
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I 1 lo 6l r ‘ s l dr = * n T n(s) , (n>l) , |s|<l 


(64) 


to regularize the system. From (62), (61) and the orthogonality of T n (r) it may be shown 
that Cj o =0, (j=l,2,3). The functional equations resulting from (60) and (62)-(64) may be 
solved by “collocation”. Thus, by specifying s (60) and (61) may be reduced to the following 
algebraic system: 


N J 


E £ c jn H Un (s im ) = gj(s im ) , (i=l,2,3; m=l,2,...,Nj) . 
J=1 n=l 


(65) 


In (65) rapidly converging results are obtained by selecting s jm as the roots of the Chebyshev 
polynomials Tj^ (s) as follows: 


s im = cos 



(i=l,2,3; 


m=l,...,Nj) . 


( 66 ) 


After solving the integral equations, the stress intensity factors defined by (46) and (47) may 
be obtained as 


k 3 ( a j) =^i2 a . /ij(x)j2(x- aj ) ^j(x) = /ij(aj)^-^- J Fj(-l) , 


b;-a: 


k 3 ( b j) = ^j(x)J 2(bj-x) ^j(x) = F,(l), (j=l,2,3) . 


(67a, b) 


Also, from (21), (59) annd (62) the crack surface displacements are found to be 
Wj(x,+0) = [ ^(x)dx = - ^ E B c in U n . 1 (s)-Jl-s^ , 

J a i 1 


= ^ ^ ■ o=w> • 


( 68 ) 


The results given by (59)-(68) are valid for all disconnected crack geometries, including 
the cracks terminating at “interfaces” x=0 and x=h. For cracks crossing the interfaces, since 
the intersection point of the crack and the interface is not a singular point, the problem can be 
solved by treating the connected crack as a single crack. For example, if the medium contains 
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a single crack along y=0, a 2 <x<b 1 (34) and (27) may be replaced by 


r b l 


k(x,t)0(t)dt = 7T r(x) , a 2 <x<b, , 

(69) 

> a 2 


r b l 


<^(t)dt = 0 , 

(70) 

) a 2 



where 

<{> 2 (t) , a 2 <t<0 , 

*(t) = I (71) 

^i(t) ? 0<t<b 1 , 


r 2 (x)//i 2 , a 2 <x<0 , 

r(*) = { (72) 

^l( x )/W( x ) i 0<x<b x , 


k(x,t) = ^ + k 22 (x,t) , a 2 <x <0 , a 2 <t<0 , 
k(x,t) = k 21 (x,t) , a 2 <x<0 , 0<t<b 1 , 

k(x,t) = k 12 (x,t) , 0<x<b x , a 2 <t<0 , 

k(x,t) = + k u (x,t) , 0<x<b 1 , 0<t<b 1 . 


(73) 


The problem may then be solved by following the numerical technique described by (59)-(68). 
Specifically, by selecting <^(t)=f(r), a 2 <t<b l9 -l<r<l, as in (62), the stress intensity factors 
and the crack opening displacement may be obtained from (67) and (68). 

Some calculated results are shown in Tables 1-5 and Figures 2 and 3. From the 
geometry of the medium it is clear that an external load giving rise to only a x z in the crack 
region would have no effect on the crack problem; that is, for this case the relative crack 
surface displacements and the stress intensity factors would be zero. In the problem under 
consideration from a practical viewpoint perhaps the most important external load would, 
therefore, be a uniform shear strain jy Z =j 0 applied to the composite medium away from the 
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crack region. Thus, the perturbation problem will be solved under the following crack surface 
tractions; 

T i(x) = ■ <T yz( x >0) = -To^j(x) . a j <x<b j , (i=l,2,3) . (74) 

Table 1 shows the stress intensity factors for two different material pairs containing a 
single crack in either material 2 or material 3. In both cases one crack tip is assumed to touch 
the interface x=0 or x=h. The table shows the normalized stress intensity factor defined by 

k(dj) = k 3 (dj)/r 0 fa , r 0 = ^ 2 y Q (75) 

where 2a is the crack length and dj designates the crack tip. For a/h— ►O, as expected, the 
stress intensity factors approach the corresponding homogeneous results, namely 

k 3 (aj) -► k 3 (bj) -4 /ijTofa , (i=2,3) , (76) 

where a 2 =-2a, b 2 =0, a 3 =h, b 3 =2a. For h/a— ►O one would expect to recover the solution of 
the corresponding crack problem in a piecewise homogeneous material (see Erdogan and Cook 
(1974)). In this case near the crack tip terminating at the interface the stress state has the 
behavior <r jz ~r’ a , (i=x,y), where r is the distance from the crack tip. Since for 

the crack in material 2 and a<0.5 for the crack in material 3. In the limit as h/a— >0, 
therefore, k(0)— >00 and k(h)— >>0. The stress intensity factors at the far ends of the cracks 
shown in the table are obtained by using the solution given by Erdogan and Cook (1974). For 
p 2 =4p 3 results showing the trends for small values of h/a are also displayed in Fig. 2. 

For a crack in material 1 with the crack tip at b x =h and having a length of 2a<h or 
2a>h, the stress intensity factors are given in Table 2. The results for the similar problem 
with a crack tip at x=0 are shown in Table 3. Again, it is expected that for a/h— ►O 
k 3 (dj)— ►/ij(dj)7o'fa. Observing that k 3 is normalized with respect to /i 2 7o4*a, this is seen to 
be the trend in Tables 2 and 3. For h/a— >0 the stress intensity factors are seen to approach 
the limiting values given in Table 1 for the piecewise homogeneous material. 

In applications the results given in Table 1 would correspond to a subcritically 
propagating crack one end of which is arrested by a relatively tough interfacial zone. If the 
interfacial zone is less resistant to crack growth than the adjacent materials, then the crack is 
more likely to initiate in the interfacial zone, propagate through, become arrested by the 
tougher of the two materials, and grow into the other. Such a process may be studied by 
using the results similar to that given by Tables 2 and 3. 

The stress intensity factors for a crack of length 2a centered at x=h/2 with2a^h are 
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given in Table 4. Note that for a/h— +0, as expected, we have 


k 3 (ai) -» k 3 (b 1 ) — /i 1 (h/2)7o 1 Ta , (77) 

that is, k— >0.5 for /i 2 /fi 3 =4 and k— ►O.l for 3=100. Similarly, for a/h — *oo the stress 
intensity factors at the crack tips approach the corresponding homogeneous half plane 
solutions 4 . For a crack embedded in material 2 or material 3 some results are given in Table 
5. In order to compare the results found in this paper with the limiting results obtained by 
Erdogan and Cook (1974) and by Kaya and Erdogan (1987), in Table 5 /i 2 /^3 * s assumed to 
be 23.08 (roughly corresponding to an aluminum-epoxy pair). 

Finally, Fig. 3 shows some sample results for the crack opening displacement 
w(s)=w(x,0) where s is the normalized coordinate defined by (59) with s=^l corresponding to 
the crack tips. Three examples shown correspond to h/a=0, 0.25 and 1, 2a being the crack 
length. In each case the centerline of the interfacial zone coincides with the quarter point of 
the crack, and the boundaries of the interfacial zone are indicated by solid vertical lines. For 
h=0 as shown by Erdogan and Cook (1974) the function <f>(x) = (3/dx)w(x,0) is discontinuous 
at the interface. This may also be seen in Fig. 3. w(x), however, becomes progressively 

smoother for increasing values of h/a. 
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Table 1. Normalized mode 3 stress intensity factors 
for a crack of length 2a located in either material 2 
or in material 3 with the inner crack tip touching 
its respective interface: ^ • 100, k-kj/r.Ja, 

7 ° I s uniform shear strain at y-+®. 




M 2 /t* 

3' 4 



crack in mat. 2 

crack 

in mat . 3 

2a/h 

k(-2a) 

k(0) 

k(h) 

k(2a+h) 

0.2 

1.0070 

1.017 

.246 

.2483 

0.5 

1.0153 

1.039 

.241 

.2464 

1.0 

1.0254 

1.071 

.234 

.2441 

2.0 

1.0385 

1.125 

.223 

.2415 

3.0 

1.0471 

1.169 

.215 

.2400 

4.0 

1.0532 

1.206 

.209 

.2389 

5.0 

1.0579 

1.239 

.203 

.2382 

10.0 

1.0716 

1.364 

.186 

.2363 

20.0 

1.0831 

1.524 

.167 

.2351 

50.0 

1.0945 

1.79 

.143 

.2341 

100.0 

1.1005 

2.04 

.126 

.2338 

200.0 

1.1049 

2.34 

.110 

.2336 

2000.0 

1.112 

4. 

.06 

.2334 

“►CO 

1 . 1144 

CO 

0 

.2333 



w 

100 



crack in 

mat . 2 

crack in mat. 3 

2a/h 

k(-2a) 

k(0) 

k(h) 

k(2a+h) 

.2 

1.021 

1.053 

.00951 

.00980 

.5 

1.042 

1.120 

.00896 

.00963 

1.0 

1.065 

1.214 

.00831 

.00947 

2.0 

1.092 

1.362 

.00747 

.00932 

4.0 

1.121 

1.584 

.00651 

.00920 

10.0 

1.156 

2.021 

.00521 

.00910 

100.0 

1.221 

4.3 

.00256 

.00903 

-►CO 

1.331 

CO 

0 

.00902 
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Table 2. Normalized mode 3 stress intensity factors 
for a crack of length 2a with the right crack tip 
touching the 1-3 interface and the left crack tip 
either in material 1 (2a/hil) or in material 2 
(2a/h>l) : “ 4 - 100 . k-kj/r.Ja, r„- 7o/* 2 > is 
uniform shear strain at y-»®. ^ 


2a/h 

k(h-2a) 

- 4 

k(h) 

V M 3~ 

k(h-2a) 

100 

k(h) 

.2 

.316 

.254 

.0216 

.0105 

.4 

.399 

.258 

.0469 

.0108 

.6 

.506 

.262 

.1034 

.0110 

.8 

.647 

.267 

.234 

.0112 

1.0 

.859 

.275 

.590 

.0116 

1.2 

.963 

.288 

.816 

.0131 

1.5 

1.028 

.310 

.966 

.0158 

2.0 

1.077 

.343 

1.092 

.0202 

3.0 

1.115 

.395 

1.202 

.0279 

5.0 

1.135 

.466 

1.282 

.0397 

10.0 

1.141 

.568 

1.336 

.0598 

20.0 

1.138 

.679 

1.361 

.0866 


Table 3. Normalized mode 3 stress intensity factors 
for a crack of length 2a with the left crack tip 
touching the 2-1 interface and the right crack tip 
either in material 1 (2a/h£l) or in material 3 
(2a/h>l) : “ 4 > 100 > k “ k 3 / r oJa, t.- 7.A* 2 ’ is 
uniform shear strain at y-*». 


2a/h 

V'V 4 

k(0) k(2a) 

V"3 _ 

k(0) 

100 

k(2a) 

.2 

.980 

.788 

.919 

.445 

.4 

.955 

.617 

.820 

.188 

.6 

.926 

.480 

.728 

.0774 

.8 

.894 

.369 

.651 

.0312 

1.0 

.859 

.275 

.590 

.01163 

1.2 

.821 

.246 

.541 

.00865 

1.5 

.773 

.233 

.486 

.00792 

2.0 

.712 

.227 

.423 

.00783 

3.0 

.635 

.225 

.347 

.00806 

5.0 

.551 

.226 

.271 

.00838 

10.0 

.460 

.229 

.194 

.00868 

20.0 

0 

.233 

0 

.00901 
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Table 4. Normalized mode 3 stress intensity factors 
for a crack of length 2a centered at h/2 in material 
1. Crack is either entirely in material 1 (2a/h<l) or 
crosses both interfaces (2a/h>l) : ^0/^3 ~ ^ > 100, 
k-k^/To-Ja, r 0*70^2* 7o is uniform shear strain at 

y-Ko , 


P 2 /^3“ 4 m 2^3“ 100 

2a/h k(h/2-a) k(h/2+a) k(h/2-a) k(h/2+a) 


.2 

.4 

.6 

.8 

1.0 

1.4 

2.0 

3.0 

5.0 

9.0 

20.0 

200.0 


.553 
.611 
.676 
.753 
.859 
.927 
.956 
.973 
.985 
.992 
.996 
.9997 
1 . 0000 


.449 

.401 

.357 

.316 

.275 

.255 

.2493 

.2477 

.2478 

.2485 

.2492 

.2499 

.2500 


.138 

.189 

.263 

.379 

.590 

.758 

.842 

.899 

.941 

.968 

.985 

.999 

1.000 


.0685 

.0454 

.0296 

.0190 

.01163 

.00938 

.00913 

.00926 

.00950 

.00970 

.00986 

.00999 

.01000 
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Table 5. Normalized mode 3 stress intensity factors 
for a crack of length 2a located in either material 2 
or in material 3. The crack orientation with respect 
to its respective interface is constant and the 
position of the other interface is varied by changing 
2a/h: m 2 /^ 3 “ 23 - 08 . k-k 3 /rJa, r.-7./* 2 > is uniform 
shear strain at y-*». J 


- (b«+a») 

* 4 If 

)0 

k(b 2 ) 

Crack in material 
1.15 

k(a 2 ) k(b 2 ) 

2.00 

k(a 2 ) k(b 2 ) 

b 

2a/h 

. V 

2* a 2 

k(a 2 ) 

1. 

1.0505 

1.154 

1.0390 

1.0869 

1.0159 

1.0244 

2. 

1.0742 

1.267 

1.0541 

1.0869 

1.0196 

1.0314 

5. 

1.1081 

1.507 

1.0717 

1.2023 

1.0228 

1.0379 

10. 

1.1324 

1.777 

1.0810 

1.2470 

1.0241 

1.0406 

20. 

1.1539 

2.143 

1.0868 

1.2791 

1.0248 

1.0422 

50. 

1.177 

2.81 

1.0909 

1.3034 

1.0253 

1.0432 

100. 

1.192 

3.5 

1.0923 

1.3125 

1.0254 

1.0435 

1000. 



1.0937 

1.3213 

1.0255 

1.0438 

-*co 

1.2558 

OO 

1.0938 

1.3223 

1.0256 

1.0438 




Crack in material I 



b_+a 

«-2h 






J 

- 1.00 

i 

15 

2 

00 

b 3* 

a 3 






2a/h 

k(a 3 ) 

k(b 3 ) 

k(a 3 ) 

k(b 3 ) 

k(a 3 ) 

k(b 3 ) 

1. 

0.0377 

.04146 

.03997 

.04182 

.04231 

.04266 

2. 

0.0346 

.04080 

.03843 

.04133 

.04203 

.04252 

5. 

0.0295 

.04009 

.03640 

. 04084 

.04178 

.04239 

10. 

0.0253 

.03975 

.03520 

.04061 

.04167 

.04235 

20. 

0.0213 

.03953 

.03439 

.04047 

.04161 

.04232 

50. 

0.0165 

.03939 

.03380 

.04039 

.04158 

.04230 

100. 

0.0134 

.03934 

.03359 

.04036 

.04156 

.04230 

1000. 


.03929 

.03338 

.04033 

.04155 

.04229 

-XX) 

0 

.03929 

.03336 

.04033 

.04155 

.04229 
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Fig. 3 Normalized crack surface displacement in bonded materials with a nonhomogeneous 
interfacial zone. In each case shown the crack goes through the zone and the 
centerline of the zone coincides with the quarter-point of the crack, s=-0.5, where s 
is the normalized coordinate defined by x=s(b 3 -a2)/2 + (b 3 +a2)/2 (Fig. 1). 



APPENDIX A 

Derivation of R 1 (a),...,R 4 ( a ) that appear in (29) 


Substituting from (12)-(18) and (28) into (7) and (8), inverting the sine transforms, 
changing the order of integrations, and evaluating the inner integrals, R lv ..,R 4 defined by (29) 
may be expressed as 


«.<•> - - a r 

J -oc 

• r°° 

R 2<“)=* 

J -OO 

= -a r 

J -cx 

. r°° 

R 4 (o) = * 

J -OO 


aA 1 (^)dg 
, a 2 +0 2 +iPO 

qgA 1 (g)dg 

a 2 +0 2 +i00 


+ A 


J -o 


- -i- f°° 

2jt I 

J -C 


aA 2 (^)d0 

a 2 +0 2 


+ ( c 2* c l) 


a0A 2 (0)d0 

a 2 +7 2 


1 


a ’ 


aA 1 (g)e~ igh dg 

i a 2 +0 2 +i/?0 

a0A 1 (0)e igh d0 

a 2 +0 2 +i/?0 


+ 2* J.„ „ 2 +* 2 + ( 3 l) “ ' 


-Af 


a0A 3 (0)e' igh d0 


oo a 2 +0 2 


(Al-4) 


We now substitute from (22) into (Al-4), change the order of integrations, and evaluate 
the inner integrals by using the following results: 

f = **<“, .< 0 , 

J -OO a I V 


A6 1 


— % — 5 - d0 = % (e at -l) , t<0 , 

0(a 2 +0 2 ) a 2V 


i0t 


^ d * 


7T -M 


a^-\-0 z -\-\f}0 


al e > t>0 > 


r 00 J0t 

e 

-oo 0(<> 2 +0 2 +i/?0) 




-11,1 


i0(t-h) 
e 

a 2 +0 2 +i/?t 


d0 


= % e" 2(h ' l) , t<h , 


28 


i0(t-h) n 2 (h-t) 

1 — s — d0 = -i*'[ e -5TnT 1- 2 ] » t<b ’ 

6(a 2 +6 7 +il36) a l n 2 a 2 


f oo W(t-y) 


a 2 +0 2 


d9 


= l ^ , »>h , 


i0(t-h) 


-a(t-h) 


- * >h 


0(a 2 +O 


(A5) 


Thus, by also using (24b, c), equations (Al-4) may easily be shown to reduce to 

R l(“) = -2^[ 6 Dl ^l (t)dt * ^ la 2 '“'MW* ’ 


R 2 (a) = - 2 ^- | e 1 <^ x (t)dt + 2 } a eat< ^2( t )^ t > 


r 3(«) = -2^r 2 £ • Mh '*w* + i , 


R 4 («) = - 5 ^ £ e n2(h '*V 1 (t)dt + 1 e- o( ‘- h) d 3 (t)d t 


(A6a-d) 
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APPENDIX B 


Expressions of the kernels ky , (i j=l,2,3) that appear in (34) 

To examine the asymptotic behavior of the kernels defined by (32) and (33) we observe 
that for |or| — *oo 

= -M\ 1+7F — l"l . <*i = “[1+(^) 2 ] 1/2 , 


n l 




A — -4a 2 . 


(Bl) 


Thus, after separating the terms that become unbounded at x=t, x=0=t, and x=h=t, the 
kernels may be expressed as follows: 


lim 

y-*+0 


Su 


hjj(x,y,t) = ti + k u( x »*) 


(B2) 


ku(x,t) = - | log|t-x| + h(x,t) + g(x,t) , 


(B3) 


QQ ^ It 

h(x,t) = [Co + logA + £ (-l) k 2i^2k)! ( t_x ) 2k J 


4- f [(l+^o ) 1 ^ 4 cos k - l]sina(t-x)da 

J o or * 

+ f (l+^o) 1 ^ 4 sin I cosa(t-x)da 

Jo <* * 

+ | [(l+^) 1/4 sin | - ^]cosa(t-x)da , 

A 


(B4) 


6 = tan' 1 (|) , 


(B5) 
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. , f(t-x) 

g(x,t) = 


/ f°° a 2 , -a 1 (2h+t-x) -a 1 (2h-t+x) 

Sj ( a - ni )( a + n 2>[-iH e 'ni e ]da 


, f A a 2 /? r -Oi(t+x) -a 1 (2h-t-x) 

+ J„ ^ [e + e 


+ 


00 r o 2 /? -o x (t+x) 
A 6 



(t+x) 


]da 


f°° -o 1 (2h-t-x) _£ e -a(2h-t-x) 

J A 4a 


+ | Ei[-A(t+x)] + | Ei[-A(2h-t-x)] l , 


(B6) 


K 12 (x,t) = I I- + f A {^[-(a+n 2 ) 
J o v 


-2aih nix . x 

e 1 e at e 1 


at e n 2 


x ]+e -a (x - ,) . 


da 




■ 2 a 1 h /vt Hi x . . /vt 1 

e 1 e at e 1 +(a+n x )e at e 


l 2 x i 1 -«( x -t) P -a(x-t)i 
J+e *4a e J 


da 


P 


- J Ei[-A(x-t)] , 


(B7) 


K 13( x >t) = {?f K«-“2> 


n*,(h-x) -a(t-h) , . -2aih -a(t-h) n 2 (h-x) 

ie 1V e v ; +(a-n 1 )e 1 e v J ] 


- e 


-a(t-x) 


}da + J~{f [-(a-n 2 > 
A 


-n 1 (h-x) -a(t-h) 
e e 


-2a x h -a(t-h) -n 2 (h-x) 1 ^-a(t-x) -a(t-x)'! ^ 

’4a e J 


, x *4Ui it -m t-u J 'lin l J. - 

+ (a-n 1 )e 1 e v J e ']-e 


- | Ei[-A(t-x)] , 


(B8) 
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K 21 (x,t) = ^+^{^[-^V ni \“ + ^V“' h e^( h - t )«“ ] .e-“ ( ‘- X) }d a 
+ J°° pa 2 [ (gHhnj^ -ni^ttx , (a+n^) ^^ n 2 (h-t) 

e ax ]e -a ( t- x ) 


- ^ e X ^}d« * f Ei[-A(t-x)] , 


(B9) 


rs , .v f\ af3 *2a!h o(t+x) 
K 22 ( x ,t) = [--£ (1-e )]e 'da 

J o 

+ f[-f (l-e- 2 “‘ h ) ■ ^]e a(t+1 ‘ ) d« 

* A 


- | Ei[-A(t+x)] 


(BIO) 


K 31 (x,t) = •^+| 0 °|2g? [ (^iJ e n 2^ e " n l t e ' o (x'h)_(^n 2 ) e ii 2 (h-t)^-a(x-h)j 


+ e 


a(x-t) 


}d„ + |”{2g?[(^i) e ”2 h e -"iV“< x 


-h) 


(a-n 2 ) n 2 (h-t) -a(x-h) -a(x-t) 0 -a(x-th . 

~ n 2 e e J + e * 4^ e / d « 


- | Ei[-A(x-t)] , 


(BH) 
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K 33 ( x »t) = f [-% (1-e 1 )] 

j o 


-2a x h -a(t+x-2h) 


da 


+ 


f[-f (l.e- 2 “‘ h ,.fje-“(‘ + - 2h )d« 

* A 


- | Ei[-A(t+x-2h)] 


(B12) 


In equations (B4-B12), C o =0.5772156649 is the Euler’s constant, A is any positive constant 
and the Exponential Integral Ei[-Az] is defined by (A>0, z>0) 

f°° 

Ei[-Az] = - \ e" az da 

J A 


= C 0 + logA + logz + 


f Az .-;.i 

Jo 


da 


oo 

= Co + logA + logz + Y 

k=l 


(-l) k A k z k 

k(k!) 


(B13) 


where rapid convergence is insured by selecting 0<A<1. 
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